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Abstract 

Polydisperse mixtures are those in which components with a whole range of sizes are present. It is shown that the 
fluid phase of polydisperse hard spheres is thermodynamically unstable unless the density of large spheres decreases 
at least exponentially as their size increases. The instability is with respect to the large spheres crystallising out into 
multiple solid phases. 



PACS: 64.10.+h, 64.75. +g, 82.70.Kj 

Mixtures of hard spheres in which spheres with a wide 
range of diameters are present are a good first model of 
emulsions. Emulsions are suspensions of droplets of oil or 
fat in water; milk is perhaps the most familiar example. 
The droplets of an emulsion interact via a short ranged 
repulsion, which is well represented by a hard-sphere in- 
teraction. They are typically present with a wide range 
of diameters: from 0.1 to a few micrometers Mix- 
tures in which a continuous range of sizes are present are 
termed polydisperse [0. They are much less well under- 
stood than systems which contain only one or two compo- 
nents. For example, the phase behaviour of single com- 
ponent hard spheres JB| has been understood for thirty 
years: the fluid phase is stable up the point where the 
spheres occupy a little less than half the volume of the sus- 
pension, there is then a first order transition to a solid. 
In contrast there are no phase diagrams known for poly- 
disperse hard spheres. Below, we examine polydisperse 
spheres with particular emphasis on the largest spheres. 
We show that unless their density decreases at least ex- 
ponentially with increasing size, they crystallise out of the 
mixture at all densities. The mixture is then never sta- 
ble as a single fluid phase. The crystallisation is driven 
by a depletion attraction ||[7| between the large spheres, 
due to the smaller spheres. Depletion-induced separation 
of the largest spheres has been observed in emulsions || 
but there the floating of the droplets to the surface due 
to gravity complicates the situation. Our demonstration 
applies to spheres at equilibrium. 

Specifying a polydisperse mixture requires specifying 
the number density of spheres of every size. This is done 
with a distribution function x(a) Q. The number den- 
sity of spheres with diameter a is then px(a)da, where p 



is the total number density of spheres. Although our final 
result will apply to a whole class of distribution functions 
we choose a specific function for definiteness and because 
it is widely used to describe emulsions pl^| and powders 
||. The distribution is called the log- normal distribution, 
and it is defined by 
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where a is the mean diameter and w is the standard devi- 
ation in units of a. Note that for this distribution there is 
no upper limit on w; its lower limit is zero and corresponds 
to a one component system. In the canonical ensemble a 
polydisperse mixture of hard spheres is completely speci- 
fied by x{a) and the total number density of spheres p; for 
hard spheres the temperature is not a relevant variable. 

In order to make progress in understanding a polydis- 
perse mixture of hard spheres with a broad distribution 
of sizes we distinguish between the spheres with diameters 
close to or less than the average diameter a and spheres 
with much greater diameters. Due to the large differences 
in size and number density of these two sets of spheres we 
will treat them differently. Not only is the number den- 
sity of the large spheres much less than that of the spheres 
with diameters near a but the fraction of the fluid's volume 
they occupy is much less. This is so because x{a) 7 Eq. (|l|), 
decays much faster than cr~ 3 at large a. 

Consider the very large spheres of the distribution, 
those with diameters a 3> a. These spheres are immersed 
in a 'sea' of spheres much smaller than themselves, for each 
large sphere there are many spheres with diameters of the 
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same order as a or smaller. These smaller spheres induce 
an effective attraction between the large spheres of the 
polydisperse mixture: the well-known depletion attraction 
first described by Asakura and Oosawa || . This effect has 
been extensively studied both theoretically R |9|-|T^1 and in 
experiments on colloids which accurately model (polydis- 
perse) hard spheres 

The depletion attraction is entropic in origin (it cannot 
have any other origin as in hard spheres there is no energy 
of interaction and so there is nothing but entropy). When 
two spheres approach each other then the volumes they ex- 
clude to the other spheres overlap. Thus the volume this 
pair of spheres denies to the other spheres decreases and 
so the volume available to the other spheres increases, in- 
creasing their entropy. This is particularly pronounced for 
a pair of large spheres surrounded by many small spheres, 
then when the large spheres touch the entropy of very many 
small spheres increases. 

For a pair of spheres of diameter a immersed in an ideal 
gas of spheres all of diameter a' <C a the range of the de- 
piction attraction is a + a' . The attraction increases from 
zero when the surfaces of the large spheres are a' apart to 
a maximum when the large spheres touch. The strength 
of the effective attraction can be measured by its value at 
contact divided by the thermal energy kT, u. This is the 
increase in the entropy of the small spheres of size a' when 
a widely separated pair of spheres of size a is brought into 
contact with each other. It is given by H[t| 

li = -p s v ov , (2) 

where p s is the density of the small spheres. Each large 
sphere excludes the smaller spheres from a spherical vol- 
ume of diameter a + a' , shown in Fig. 1 by the thick lines 
around the spheres. When two large spheres are touching, 
the two volumes which they exclude to the small spheres 
overlap. The volume of overlap of these two volumes is 
v ov . In Fig. 1 we see that this volume is equal to that of 
two caps, each an end of a sphere of diameter a + a' and 
of height a' /2. We are considering the limit of small a' ju 
and so the caps are very flat. Then the height of one of the 
caps a distance x from a line drawn between the centres of 
the two large spheres is (<x'/2)(l — x 2 /r 2 ) where r is the 
radius of a cap at its base; r 2 — aa' /2. The total volume 
of the two caps 
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Now, for polydisperse spheres distributed according to 
Eq. (|IJ) the number density of spheres with diameter a' is 
px(a')da' and so for a pair of spheres of diameter a a c 
the depletion attraction due to spheres with diameters < a c 
is 

u(a) = —~rpa / x(a')a' da' , (4) 
4 Jo 

which depends on the cutoff a c . However, if a is suffi- 
ciently large that a <C a c <C a then the integral depends 
only weakly on a c because x(a')a' 2 is small for values of 



a' > a c . Indeed we can extend the upper limit of inte- 
gration to infinity without introducing a significant error. 
The fact that we can do so justifies our splitting of the 
distribution into two parts. Then we have 
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where we have used the relation 77 = (ir/6)pcf 3 (l + w 2 ) 3 
which holds when x(a) is given by Eq. (|lj). The physi- 
cal content of this approximation is that very large spheres 
only notice spheres with diameters around a and less; the 
density of the larger spheres is too small to add significantly 
to the depletion effect. Note that the attraction increases 
linearly with the size of the spheres a. 

By using the idea of a depletion attraction we have re- 
duced our polydisperse mixture to the large a tail of the 
distribution interacting via an effective interaction which is 
the sum of a hard-core interaction plus the short ranged at- 
traction of Eq. (||) . The attraction of Eq. (||) favours con- 
densed phases where the large spheres are within the range 
of the attraction of each other. Competing against this at- 
traction is the translational entropy of the large spheres, 
which favours dilute phases. This competition is the same 
as that involved in the vapour-liquid transition of a simple 
substance such as argon. The translational entropy of the 
large spheres is just that of an ideal gas mixture, so per 
large sphere of size a it is B 
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where the second expression is obtained by substituting Eq. 
(|l|) in the first and the constant is a function only of w. The 
entropy Sp increases with sphere size a because the den- 
sity decreases. But it only increases as the square of a log, 
which is a slower than linear increase. 

The stability of the dilute fluid phase with respect to 
condensation into a phase in which the density of large 
spheres is much higher, is determined by the relative en- 
tropy of the dilute and condensed phases. Therefore, we 
require the entropy of the condensed phase. The range of 
the depletion attraction is ~ a '. For spheres of diameter 
a ^> (7, this is very small in comparison to the size of the 
sphere, a. In the a — > 00 limit, the ratio of the range of 
the attraction to the size of the sphere, a /a, tends to zero. 
This is the sticky-sphere limit introduced by Baxter ]17j . 
Thus the large spheres are all near this sticky-sphere limit 
and as a diverges the spheres tend towards the limit. Stell 
has shown [jl8],[l!| that as the strength of the attractions 
is increased, a fluid of sticky spheres does not condense to 
form a liquid but collapses to form a close-packed solid, see 
also Refs. p0[|21{ ]. So, we look not for condensing of the 
large spheres into a dense, liquid-like, phase but for collapse 
into a dense solid. We therefore require the entropy in the 
dense solid phase. By dense we mean sufficiently close to 
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the close packed density that the sphere is within the range 
of the depletion attraction of its neighbours. 

The entropy per large sphere, of a dense solid phase 
has two parts. The first is the entropy associated with the 
motion of the large sphere, and the second is the entropy 
gain of the small spheres when a large sphere is brought 
close to twelve neighbouring spheres, as it is in a dense 
face-centred-cubic or hexagonal-close-packed lattice. The 
first part is easily obtained from a cell theory [^j],^]. This 
assumes that each sphere is restricted to a cage formed by 
its neighbours, which are taken to be fixed at their lattice 
positions. For a solid with a lattice constant a, the centre 
of mass of a sphere can move a distance ~ (a — a) from 
its lattice position without bumping into any of its neigh- 
bours. The solid is formed due to attractive interactions so 
the spheres must be close enough to each other to attract 
each other throughout the cage which its neighbours form. 
For this to be true a must satisfy a < a + c?/2, and so 
a — a = ca, where c is a parameter in the range (0, 1/2). 
For monodisperse spheres the entropy is simply the loga- 
rithm of the volume (ccf) 3 , which is available to the centre 
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When the spheres are polydisperse the situation is more 
complicated. There is an upper limit to the range of sizes 
of spheres a single solid phase can tolerate p3|-p^|. A lat- 
tice can only accommodate spheres up to about its lattice 
constant a in diameter; larger spheres cannot fit into the 
lattice position without overlapping with their neighbours. 
Spheres with diameters less than ~ (a— a) are so small that 
they cannot be within the range of the depletion attrac- 
tion of all their neighbours p]J . This means that the large 
spheres cannot all crystallise into a single solid phase. In or- 
der to crystallise they first fractionate into many fractions, 
each containing spheres with only a narrow range of diame- 
ters. The fractions can then crystallise individually to pro- 
duce separate solid phases, each containing spheres of a dif- 
ferent size. The combined fractionation-and-crystallisation 
of polydisperse sticky spheres is discussed in Ref. [El] . The 
range of diameters is roughly a — a = ccf. This width of 
distribution contributes an amount ~ ln(c<7) to the entropy 
of the solid |,|2|. 

The parameter c will be determined by a competition 
between the depiction attraction tending to reduce it and 
the motion of the large sphere which tends to increase it. 
However, our results are not sensitive to the exact value of 
c and so we merely take it to be much less than one. Then 
the depletion attraction is almost equal to its value at con- 
tact, Eq. (||), and the gain in entropy of the small spheres 
per large sphere which solidifies is closely equal to minus 
six times Eq. (|5|). Then, the entropy per large sphere of a 
solid phase of large spheres of size a 

sk{o) ^ 41n(cCT) - 6u(a) er » ct, c < 1. (8) 

The solid phases are much denser than the fluid phase 
and so can be formed without increasing the volume occu- 
pied by the system. Thus if the solid phases have a higher 
entropy than the dilute fluid, then the dilute fluid cannot 



be the equilibrium phase, as the entropy can be increased 
at fixed volume by forming the solid phases. So, we now 
compare the entropies of the fluid and solid phases. The 
entropy change per large sphere As, when spheres of av- 
erage diameter a separate out from an ideal gas to form a 
crystal phase with a polydispersity of order ccf is Eq. (||) 
minus Eq. (Jfy 
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When the volume fraction r\ is non-zero, this is positive for 
sufficiently large a. In fact it is positive for any x(a) which 
decreases more slowly than exponentially with a. There- 
fore, any fluid phase of hard spheres with a distribution 
which decreases more slowly than exponentially is unstable 
with respect to the largest spheres crystallising into solid 
phases with narrow polydispersities. A single solid phase 
can contain only a narrow slice, of width a fraction of a, of 
the original distribution x(a) but spheres with diameters 
ranging from infinity down to some large but finite limit 
crystallise. Thus, an infinite number of solid phases form; 
each phase with a different range of sphere sizes. For the 
sake of clarity, when we say that the fluid phase is unstable 
we mean that the solid phases have higher entropies. The 
fluid phase will however be metastable, i.e., stable with re- 
spect to infinitesimal perturbations. 

Very recently, Cuesta |27j has shown that within the 
Boublik-Mansoori-Carnahan-Starling-Leland (BMCSL) Q 
approximation polydisperse hard spheres with a log-normal 
distribution with a sufficiently large standard deviation w 
have a spinodal. Warren has also found a spinodal within 
the BMCSL approximation pS) ]. A spinodal is where the 
fluid phase becomes unstable with respect to an infinites- 
imal density and composition fluctuation. The difference 
between Cuesta's result and ours is probably due to one or 
both of two factors. The first factor is the nature of the 
transition we have found. It is very strongly first order and 
so the transition occurs much before the spinodal. The sec- 
ond factor is the accuracy of the BMCSL approximation. 
It may be poor when there are spheres of widely different 
sizes present [f29f . 

In comparing our result with experiment it should be 
remembered that in an emulsion there will be some upper 
size limit, beyond which there are essentially no particles. 
Obviously, the number of phases which separate out is then 
not infinite. In addition, at sufficiently low volume frac- 
tions the fluid phase of the emulsion will be stable. The 
fluid phase can be destabilised by adding small spheres to 
the distribution, so increasing the strength of the deple- 
tion attraction. For emulsions, micelles can be added and 
indeed this is done in Bibctte's || procedure for fraction- 
ating emulsions. 

To summarise, polydisperse hard spheres with a non- 
zero volume fraction and distributed according to a dis- 
tribution function which decays more slowly than expo- 
nentially are thermodynamically unstable. Spheres above 
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some lower size limit crystallise due to the depletion attrac- 
tion induced between them by the presence of the smaller 
spheres of the distribution. We have not determined this 
lower limit but it is much larger than the average size a. 
The solid phase can only tolerate a very limited polydisper- 
sity [pl],[23|,p5|,pq| and the range of spheres which crystallise 
is from this lower limit to infinity. Thus, the number of solid 
phases which form is infinite. This seems surprising at first 
but in the large a tail of the distribution the depletion at- 
traction is increasing more rapidly than the translational 
entropy in the fluid phase and so there is no upper limit 
to the sizes of spheres which crystallise. The lower limit 
to the sizes of spheres which crystallise is clearly finite as 
a solid phase will form so long as it reduces the entropy 
by any non-zero amount. The sub- linear increase with a 
of the translational entropy in the fluid is the crucial fac- 
tor in destabilising the fluid phase, ft inevitably leads to 
the fluid being unstable when the attractions grow linearly 
with sphere diameter. 

Finally, we conjecture that the instability we have found 
is not restricted to spheres or to attractions which arise 
from depletion. Consider a general polydisperse fluid with 
a number density px(l)dl of elements of size I then the 
translational entropy increases as minus the logarithm of 
x(l). If the attractive energy (over kT) u(l) between ele- 
ments increases faster than lnx(l) 

dKQj > d(lns(0) 
dl dl ' 

then the attractive energy is much larger than the trans- 
lational entropy for sufficiently large I. When this is true 
we expect the mixture to be unstable with respect to the 
largest elements condensing out to form a dense phase in 
order to minimise the energy. 

It is with pleasure that I acknowledge discussions with 
J. Cuesta. 
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Figure caption 

Fig. 1. A schematic of two large touching spheres of di- 
ameter a, the shaded discs, with the volumes they exclude 
to smaller spheres of diameter a' . These volumes are out- 
lined by the heavy circles and they overlap when the large 
spheres touch. 
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